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I. INTRODUCTION 


The goal of this project is to develop a direct sequence spread spectrum binary phase 
shift keyed (BPSK) communications subsystem for the Petite Amateur Navy Satellite 
(PANSAT). Previous work on this topic primarily focused on modulator and demodulator 
(MODEM) designs that utilized various modulation techniques, such as BPSK, and four 
frequency shift keying (4-FSK). These MODEMs were designed to operate at 1200 bps, 


have a pseudo noise sequence length of 127, and have a chip rate (f,,..) of 152.4 kHz. 


hip 
With a chip rate of 152.4 kHz, the null to null bandwidth of the transmitted signal was only 


2b nip 


= 304.8 kHz. The link analysis, which determined these and other parameters such as 
a required transmitter power of 7 watts, was found to be in error. Fortunately, the 
conservative link analysis that 1s included as part of this thesis provides communication 
system performance parameters that far exceed those previously discussed. However, the 
new parameters called for more bandwidth than was previously understood to be available 
from the Federal Communications Commission (FCC) and the American Radio Relay 
League (ARRL). As it turned out, the project only requested 1 MHz of bandwidth, and that 
is all the ARRL awarded, 1 MHz of bandwidth centered at 437.25 MHz. After a review of 
all documentation and applicable rules manuals, this too was found to be in error. The 
rules actually allow space operations anywhere from 435-438 MHz. So to accommodate 
the wider band signal the center frequency was moved to 436.5 MHz. 

Part of the problem was that the communications system was not being looked at, or 
designed as a system. Previous efforts, in addition to being based on faulty data, ignored 
the rest of the communications problem by concentrating solely on the MODEM designs. 
Following the discovery of the errors previously mentioned, it was recognized that the 
communications problem needed to be approached as a system. So this thesis project was 


expanded from the design of a MODEM to the design of a complete communications 


subsystem. The intent was not only to design, but to design, build, and test the complete 
communications subsystem. 

This thesis begins with an explanation and results of the most recent link analysis. 
Followed by a functional description of the communications subsystem that evolved 
during the course of this thesis work. The details of this work, such as the schematics and 
board layouts have been included as appendices. Additional information on pseudo noise 


sequences, the link analysis, and filter design may also be found in the appendix. 





II. LINK ANALYSIS 


Before attempting a design of either the transmitter or receiver, a link analysis was 
performed for both the up-link and down-link. Past experience has shown that even though 
the link budget equation is very simple, minor errors are easily made. That is why a general 
link analysis program was written using MATLAB® and included as Appendix B. Most 
of the equations used in the program came from [Ref. 1:p. 390-396]. Some of these 
equations are repeated here for clarity. This Matlab® program allows the user to quickly 
vary the system parameters and plot the results over a wide range of transmitter power. The 
results are output in both graphical and tabular formats. The tabular format provides, 
among other things, an estimate of the received carrier power at the base of the antenna. 
This estimate is then used to compute the power throughout the design. 

As one might expect, the receiver will have to process a very weak signal; just how 
weak is what needed to be determined. The answer was provided by the link analysis 
program. A main concern was that the noise added by the receiver components might 
obscure this weak signal. Fortunately, the receiver’s noise figure may be controlled though 
the proper selection of receiver components such as the low noise amplifier (LNA). Once 
the receivers noise figure is known, the carrier-to-noise-ratio may be computed. Because 
each amplifier stage amplifies the signal and noise equally, the carrier to noise ratio at the 


front of the receiver is the same as the carrier to noise ratio at the detector. 


§)- @,-@, en 


Since this ratio remains constant throughout the receiver, the noise power may be ignored 


and the signal power used exclusively in power calculations. 


Table 2-1: System parameters used in the link analysis. 


Satellite Antenna Gain Go atellite = 0 dB 
Antenna Noise Temperature T, = 280° Kelvin 





For the purpose of comparison, the transmitter power was varied from 0.2 watts to 2.0 
watts and the resulting effective isotropic radiated power (P,j,p) Computed using the 
following equation. 


Prirp = PG, (2.2) 


Parameters such as P, and G, are shown in Table 2-1. Once the Pear is known the 


received carrier power may be computed using equations (2.3) and (2.4). Equation (2.4) is 
used to determine the amount of loss due to free space. 


Cry = (Peipp) + Gap Legs (2.3) 





4 
(a 20108) ad (2.4) 


The noise component of (=) is computed using equations (2.5) through (2.8). These 
Rx 


equations are used to compute the system noise temperature. Looking at equation (2.8) it | 


is easy to see why the low noise amplifier (LNA) is the single largest 





contributor to the receivers noise figure and noise temperature and why its selection in 


terms of gain and noise figure 1s so important. 





Ney = kT ,y:Byp (2.5) 

re tat (2.6) 

T, = T,(F,-1) (2.7) 

F +. Z ie ee eee (2.8) 
° Gina GinaGy GinaGy---G, 


Since this link analysis considers only the worst case, the slant range (d) was only 


computed for an elevation angle of E = 5 using the following equation. 








a= (R+n)? +R? 2K (RH) sin( EX 4 asin = Jeos (E }} (2.9) 
180 (R+h) 180 


The results are then converted to dB’s, thus reducing the problem to simple addition. 


(Peirp) apy = log (Pripp) (210) 
Kyy = MWlog (1.38x10) = -228.6,, (2.11) 
Bee (lon (Ba) (2.12) 
4 ea) 
ee a 2 =44s eee | Se le 2.13 
& a (Peirp) apy ~ (Les) ap + ta dB~ YaB (2.13) 
E 
f=) ~ | eR (2.14) 
NE dB N dB 


E ° 
Now that is known, the probability of a bit error may be calculated using the error 


oO 


correction function as shown in the following equation. 


PL = {1-ene{ | (2.15) 


A graph of the bit energy-to-noise-ratio as a function of the information bit rate and 
transmission power is shown below. 


Bit Energy to Noise Ratio vs Information Bit Rate 


Eb/No (dB) 





10000 12000 


0 2000 4000 6000 8000 


Information Bit Rate (ops) 


Figure 2-1 Bit energy to noise ratio vs. the information 
bit rate over a power range from 0.2 to 2.0 watts. The 
lowest curve is for 0.2 watts. 


Graphs showing the probability of error are not included here, because the probability of 
error is so low, but may be found in Appendix B. Appendix B also contains the 


MATLAB® program and a set of tabular results for the data in Table 2-1. 





Hl. FUNCTIONAL DESCRIPTION OF THE COMMUNICATIONS 
SUBSYSTEM 

The following discussion is conducted using simplified block diagrams of actual 
boards designed during the course of this thesis work. Following a brief explanation of the 
overall system, each board will be discussed individually. If more detail is desired, it may 
be found in the schematics and layout diagrams included in Appendix D. If the reader 
requires areview of spread spectrum communications or wishes to see mathematical proofs 
of concepts applied here, References 5, 7, and I1 are suggested. 

Equations used here may be written differently than those found in most text books. 


In most communications text books a signal is described mathematically as 


St) = J2P cos (w,t+ 6 (t-1) — >) (3.1) 


where J2P is the power of the signal, @(t-—t) 1s the phase modulation with t representing 
a phase shift due to the channel, and 9 is the carrier phase shift due to the channel. Since 
this thesis is a practical application of communications theory and mathematics, signals 
will be described in the following way 


S(t) = Pup, COS (Wt + 6 (t- T) — >) (3.2) 


It seems more practical to write signal equations in this manner, because component 
catalogs express gain, loss, and power in units of dBs and dBms. Now gains and losses 
may be accounted for with simple addition and subtraction. 

S(t) = (Pyp,,+ Gap — Ly) cos (mt + 8 (t—t) — 9) (3.3) 
With signal power expressed in this way, the reader may easily apply these equations to 
the circuits in Appendix D, and use a spectrum analyzer to measure signal levels at critical 


points in the design. Measurements may be made by simply tuning to the carrier 


Q) e 
frequency f, = * , adjusting the span to => 2(data rate), and moving the marker to the peak 
Tt 


. . ‘ A : Z : 
of the signal waveform, which in this case is a SINC’ function. 


A. OVERVIEW OF THE COMMUNICATIONS SUBSYSTEM 
Every spread spectrum communications system must perform the same basic | 
functions. The transmitter must spread the data spectrally, up-convert to the carrier 
frequency, amplify the signal and radiate it from an antenna. Every spread spectrum 
receiver must down-convert from the carrier frequency, provide doppler compensation if _ 
necessary, amplify the received signals, detect the desired signal, track it, de-spread it, and | 
demodulate it. What could be simpler? Well, the devil is in the details. 
This communications subsystem affords the user flexibility by providing the ability | 
to: 
¢ Select between redundant transmitters and receivers 
¢ Select a desired transmitter power from a range of power levels 


¢ Evaluate numerous PN phase search strategies, and use the one that provides the 
quickest acquisition time possible 


¢ Switch between straight BPSK and spread spectrum BPSK 

¢ Compensate for doppler shift using the ground station only 

¢ Develop and test a faster acquisition circuit in the future 
Starting with a link analysis this communications subsystem was designed as a system. 
The design is modular, with each board performing a different and unique function. | 
Circuit design and component selection was based on power levels that were determined | 
from the link analysis and used to calculate the signal power throughout the design. © 


Because of the interdependency on signal power levels between each board, an iterative | 





approach to the design was employed. The design evolved, as the designer became more | 


aware of the practical aspects of circuit design and component selection. 
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Figure 3-1 is a functional diagram of the present design, where: 


¢ Board 1: provides up-conversion, down-conversion, and amplification. 


¢ Board 2: provides receiver IF amplification, automatic gain control (AGC), and 
coarse PN phase detection. 


¢ Board 3: provides fine PN phase alignment and tracking. 


¢ Board 4: de-spreads the signal, locks on to the phase of the suppressed carrier, and 
demodulates the data. 


¢ Board 5: provides modulation and amplification of PN sequences developed on 
board 6 and applied to boards 2-4. 


¢ Board 6: provides PN phase generation, search, and control logic, as well as any other 
digital circuitry used in the design. 


Though not yet complete, this design makes significant strides in the right direction. With 
its flexibility and modularity, the opportunity is there for future designers to continue its 
evolution or use it as a test bed for new and improved designs. 
B. BOARD 1: TRANSMITTER AND RECEIVER FRONT END 

Board | groups all the components with SMA connectors together. Connections via — 
SMA and 50Q2 cable eliminate printed circuit board (PCB) design problems that would | 
likely occur at frequencies above the first IF frequency of 32.1 MHz. Board 1 is designed 
to serve in either the satellite, using a fixed 404.4 MHz oscillator or in the ground station, | 
using a variable 404.4 MHz oscillator. The variable oscillator’s control input will come | 
from a program called the Kansas City Tracker that predicts the orbital track of the satellite. | 
This method of doppler compensation eliminates the need for tunable filters in the satellite | 
and places the burden of doppler compensation solely on the ground station. Another more | 
obvious way to achieve doppler compensation is to simply widen the filter bandwidths. | 
This method was ruled out because the maximum doppler shift is projected to be + 10 kHz, | 
while the information bandwidth is < 20 kHz null to null, making the noise power | 


bandwidth twice as large as the signal bandwidth. 
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The single pole four throw (SP4T) PIN diode switch allows the user to switch 
between two receivers and two transmitters, using two TTL control lines. The cavity 
bandpass filter (BPF) following the SP4T PIN diode switch serves as a preselection filter 
for the wideband low noise amplifier (LNA). This preselection filter should eliminate 
strong out of band carriers by reducing the noise power bandwidth. The LNA will provide 
needed amplification prior to mixing for downconversion to the first IF frequency of 32.1 
MHz. The inductor and capacitor (LC) BPF is then used to select the difference frequency 
from sum and difference frequencies that result from mixing. The output of the LC BPF is 
then feed to the receiver IF strip on board 2. 

The upconverter translates the 32.1 MHz IF frequency from board 5 up to the RF 
frequency of 436.5 MHz. Here, the BPF following the mixer is used to select the sum of 
the local oscillator and IF frequencies. A cavity filter is used because of the high Q and 


high center frequency required by the application. The Q 1s calculated as: 





f 436.5MH 
Q = a Ts (3.4) 
BW 2.52MHz 
Unfortunately, the size of these filters is very large when compared to other components. 
However, they are passive and have SMA connectors, as do all the components on this 


board, so they may be placed anywhere in the satellite. 


Because there is still some uncertainty over the amount of power required to complete | 


E 
the link with an acceptable - an electronic step attenuator was built into the design. | 


Detrimental channel effects such as multipath fading, are the primary concern. The 
electronic step attenuator will allow the controlling ground station to select the most 
appropriate transmitter power between 0.2 watts and 2.0 watts. The two amplifiers 


following the electronic step attenuator are expected to provide amplification up to a 
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maximum of 2 watts. However, since the project is only in the bread board stage a lower 
power amplifier was used. 

The received signal will arrive with a random PN sequence phase and random carrier 
phase. Each is indicated in the equations that follow, as t for the PN sequence phase and 
¢ for the carrier phase. For simplicity, additive channel noise is neglected. 

Oe Cos ato, (tT) 0) (3.5) 
As mentioned above, LO] may be either fixed or variable depending on the application. 


Only the fixed oscillator is considered here. 


Sto1 (0) = (Pap ~ 3a 7 Lip) ©O8 (OL9)0 (3.6) 
Sroi() = PrypmfS (@) 90 (3.7) 


After down-conversion to the first IF frequency, the following is delivered to the receiver 


IF strip on board 2. 


Sy WU) Sle rary Gos Ul Carre aa (9) ot) (3.8) 
PE) ECAB Ti), Ee LNA cond = Le (3.9) 
S(O = Pargnny (cos@ipjt+84(t- 2) = 61 (3.10) 


The upconverter for the transmitter occupies the bottom half of the board. Its input 


comes from board 5. 


S(t) = PE apm COS (Oyp)t + By (t)) (3.11) 


Naturally, the same local oscillator (LO1) is used for the upconversion to the RF 


frequency. 
SG = EaanCes (Ona b) (3.12) 
Following upconversion and amplification, the signal delivered to the antenna Is 


S,(t) = P COS (W, pt + 6, (t)) (3.13) 


t (dBm) 


P ere ya Gh (3.14) 


t(dBm) = PE dBm) = Oa Attn 


13 


With only a few minor changes, this design may operate in either spread spectrum 
BPSK and straight BPSK. The addition of a PIN diode switch and narrowband filter, in | 
parallel with the wideband filter following the mixer in the receive path, is necessary to | 


reduce the noise power delivered to boards 2 and 4. The PIN diode switch would be used | 





to select the correct filter, based on the desired mode of operation. Next, the output from 
board 6 could be switched high, to lock the bi-phase modulators on board 5 in the zero 
phase state. The controller would then only look to boards 2 and 4 for an indication of 
detection and lock. On the transmit side, no change to the filter bandwidths is necessary. 


Naturally, since the signal is to remain narrowband, the XOR gate on board 6 would be | 





switched out. So, only boards 1 and 6 are affected by this change in the mode of operation. | 


C. BOARD 2: RECEIVER IF AMPLIFICATION, AGC, AND DETECTION 





The receiver IF strip and detection circuits were placed on the same board for 
convenience because they are both relatively small circuits. However, to allow for testing | 
of a different detection circuit in the future, the two are connected via a 50Q cable. 

The receiver IF strip provides amplification and automatic gain control (AGC) prior 
to distribution to the detection, tracking, and demodulation circuits. An AGC circuit was | 
built into the design to compensate for anticipated signal fading due to multipath and the | 
satellite’s tumbling. A predetection feedback path was used, as a result of information read — 


in [Ref. 4:p. 1624-1628] and the fact that it was easier to build. 
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The wideband signal coming into the receiver’s IF strip may be expressed as 
Sa (t) = Py capmy OS (@p,t+ 9, (t-7) - 9] (3.15) 


and the amplified signal leaving the IF strip as 


S,(t) =P t+0(t-—T) -6) (3.16) 


COS (W 


B (dBm) IE} 


where Pp (apm) 1S a result of the gains and losses seen by the signal along the IF strip. 


22dB 
~ Lit (ac) — O43 ~ Lip (3.17) 


Pp (aBm) = Pa (apm + G2 + G3 4+Gy+ Gs + Gage — 





The AGC amplifier used in the circuit and expressed in equation 2.3 provides both | 


amplification and attenuation over a range of 36 dB [Ref. 6:p. (14-18)]. Naturally, the 


amount of gain provided by the AGC amplifier is a function of the voltage output from the | 


analog level detector. This voltage (V,, 


to provide a feedback voltage that varies between zero and five volts. 


Graces EQpeGe (Sus) 


A 


Once amplified, the signal is distributed to the detection, tracking, and demodulation | 


circuits via a four-way power splitter. 


The detection circuit receives its input from the four way power splitter on the same 


board via a 50Q cable to allow for the connection and testing of a different detection circuit | 
in the future. This circuit is an implementation of a single dwell detector. The concept is | 
very simple. A search algorithm, implemented using digital components, is used to shift | 
the phase of the locally generated PN sequence. This locally generated PN sequence is then | 


bi-phase modulated onto acarrier from LO2, amplified, and input to the local oscillator port | 


of the mixer 


S02 (1) = Payp C08 [®, oot + 9, (t- 2) (3.19) }} 


When the received and locally generated PN sequences come within 1/2 of a chip, | 
the received signal will be de-spread from wideband to narrowband. The narrowband | 
BPSK signal will ‘pop-up’ in the BPF following the mixer. To limit the amount of noise j 


power bandwidth, the bandwidth of this filter should kept as narrow as possible. Since data 


) is in millivolts and must be amplified in order | 








ee 


is not recovered in the detector, intersymbol interference is not a major concern. In fact, 
the bandwidth may be as low as one times data rate, to improve performance in a high noise 
environment [Ref. 5:p. 161]. To achieve a narrow bandwidth at the second IF frequency, 


a high Q filter was required. 


f 
Q = ~ = a 2713 (3.20) 
Z 








A crystal filter was the only way to achieve a Q this high. Once de-spread and filtered the 
signal may then be detected. The threshold detector receives its signal from the crystal 
filter, then provides a TTL signal to board 6, indicating coarse phase alignment. 

The autocorrelation of the received and locally generated PN sequences is a stochastic 
process. Thus the reason the search algorithm must dwell for a period of time, generally 
the code sequence length, before a decision can be made. If c(t) 1s a PN code sequence 


then the autocorrelation of c(t) and a shifted version of itself c(t-— 1) may be expressed as 
Z Z 2 2 Z 
R(t) = E{c(t)c(t-t)) = E{c] - (E[c])” = o -p = var— (mean) (3.21) 


Or 


R.(t) = c(t)c(t—t)dt (3.22) 


Z| 
Ce, tt |Z 


N| Zz 


where the overbar denotes the ensemble average [Ref. 7:p. 351]. The autocorrelation 
function may be viewed conceptually, with respect to threshold detector, as a magnitude 


scaler to the signal energy in the bandpass filter. 


S(t) = Ss Go) i an COS [ (@,5, — ML go) t+ 8, (t-1) - 9) (3.23) 
S.(t) = R, (7) P age cos [M,.,t+ 6, (t-T) -—9] (3.24) 
R(t) = 1 when t = €N and R,(t) = ~. when t#€N (25) 
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where e is an integer. When the received and locally generated PN sequences come to 
within 1/2 of a chip of perfect alignment, the energy delivered to the threshold detector 


R_(t)P >threshold). The level of the | 


should be enough to cause a detection ( “dene 
Cc m 





threshold may be adjusted by using a variable resistor to change the sensitivity. 
D. BOARD 3: RECEIVER TRACKING CIRCUIT 

The receiver tracking circuit is an implementation of the non-coherent double dither 
loop (DDL) [Ref. 5:p. 188-189]. This scheme utilizes two channels instead of one. The | 
use of two channels and a subtraction circuit provides the error voltage with a sign (+), thus | 
indicating the direction of phase change required to increase correlation. Dithering the | 
early and late PN sequences between the two channels averages out any DC offset or | 
imbalance inherent in the two channels. The frequency at which they are dithered is 
determined by the bandwidth of the BPF’s in the arms of the tracking circuit [Ref.5:p. 175]. 


BW 5 pr 
, (3.26) 





Dither frequency = 


The error voltage that is produced is then amplified, lowpass filtered, and fed back to the | 
VCXO of the PN sequence generator on board 6. A window comparator is used to | 
provide a lock detection signal to board 6. | 

The wideband signal that is supplied to the tracking circuit board 2 may be expressed | 
as 


Sy (t) = Pa apm) COS (Ojpyt + 84 (t-T) — 6) (3.27) a 


This signal is then amplified before splitting to provide a stronger signal to each channel. 
S.(@) = (Ppapa te Geeta — Upp) cesi@ye t + Ore 0 (3.28) | 


Saas cos [@,,,t+8,(t-t) -9] (3.29) | 


c (dBm) 
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The early and late versions of the estimated PN sequence are then used to bi-phase 
modulate the sinewave produced by LO2. This signal is then amplified to +7dBm and 
applied to the analog switches used to control the dithering. The other side of these 


switches is connected to the LO port of the double balanced mixers. 
Sion (Oe= aaa COsmone On (lssAT =) | (3.30) 
The autocorrelation function may be viewed conceptually, with respect to the analog 


level detector, as a magnitude scaler to the signal energy in the bandpass filter. 


Sa(Q) = (Po capmy 7 Lconv 7 Lit ppg) Re (4) 08 [ (Op) — M192) t+ 9, (t- 7) - 9] (3.31) 


S4(t) = Py capi Re (t) C08 [Wjpot + 8, (t- 1) — 9) (3.32) | 


The voltage V, produced by the analog level detector is obviously a function of the signal | 
| 


power applied to it. Its output varies linearly with the input, as shown in [Ref. 6:p. (5- | 


11)J. 





Vv, = ey anene (3.33) | 


As stated before, the error voltage is produced by subtracting the detector outputs of each 
channel. 


e(t) = (V,-V,)g(t-nT) (3.34) | 


This error signal will vary in magnitude at the same frequency as the dither signal. 


g(t—nT). 


Dither frequency = 





(3.35) | 


dither 
Therefore the cutoff frequency of the loop filter should be much less than the dither | 


frequency. 
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J 
Loop filter . T 
dither 


BW (3.36) 


Naturally, this circuit is not used when the communications subsystem is operating in 
a straight BPSK mode. Everything may still be allowed to function, but neither the 
feedback voltage nor the loop lock detection are used. 
E. BOARD 4: CARRIER TRACKING LOOPAND DEMODULATOR 

The demodulating circuit is an implementation of a Costas loop [Ref. 7:p. 145] and 
[Ref. 1:p. 293]. Its input comes from the four-way splitter on board 2. 


Sp (0) = Pp (apm COS (@yp yt + Oy (t- 7) — 9) (3.37) 


Assuming acquisition is successful, the locally generated PN sequence is in phase with the 


received PN sequence. So t = t and @ (t—7%) 1s a punctual version of 8 (t—1). 
ST aaa COS Onst +0: (t= 7))) (3.368) 


Since the received signal is being de-spread with a punctual version of the PN sequence, 
the autocorrelation function is unity and the signal is maximum. 


S.(t) = R, (4) Po apm COS [ (@yp, — ML 92) t+ 8, (t- 7) - 9] (639) 


Deane. = Eee) + Cin L 3.,-L (3.40) 


10° “conv. ~dB IL 

Now that the phase of the PN sequences is locked, the BPSK signal must be 
demodulated. However, before demodulation can take place, the suppressed carrier must 
be tracked. A Costas loop is a phase locked loop designed to lock on to suppressed carrier 
signals. Like the tracking circuit, this scheme also utilizes two channels. But, the signals 


in these channels are in quadrature. 


Soe) aun cos (ar t=O) (3.41) 


Shope (lt) =eanecOs (opeat+ : -3) (3.42) 
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Since they are in quadrature, the multiplication of the two channels yields their phase 
difference. This phase difference is expressed as an error voltage that is amplified, filtered 


and fed back to a 10.7 MHz VCXO on board 6. 


Si(t) = (Pogpm ~ Leonv + Geic) COS [ (Mp2 — Opp) t + 6, (t-T) - 6] (3.43) 


Sg (t) = eet ci) cos) (Oye Opa) t+ 5 +0 (t—T) a) (3.44) 


Now that the correct phase has been determined for both the locally produced PN 
sequence and 10.7 MHz oscillator, the BPSK signal may now be demodulated. Coherent 
demodulation is accomplished with an integrate and dump circuit. The integrate and 
dump is implemented with an active integrator followed by a sample & hold circuit. The 
result is a NRZ signal that may be converted to digital with a comparator. The timing 
signal needed to control the integrate and dump is be derived from the PN sequence clock. 
This is possible because each data bit is an integer multiple of a PN sequence chip. The 


output of the sample & hold circuit may then be written as 


m(t-t) = Viet m= 41 (3.45) 


where 6 = o and ce? = 1. The NRZ signal m(t—t) is then input into a comparator and 
converted to a binary message string m(t—t) = (0,1). 
F. BOARD 5: MODULATION AND AMPLIFICATION 

To reduced the number of components on boards 3 and 4 (the tracking and 
demodulation circuits), the modulation and amplification circuits of these boards were 
placed on a new board named board 5. This board later became the modulation and 


amplification board for both the transmitter and receiver. 
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The circuit that produces S, 5, and S, 55, will allow testing of a faster acquisition 


circuit in the future. 


S102(0 = (Prapm ~ Sap — Lat ~ 3ap — Lyi) €08 (@z oat) (3.46) 
Sto27a(O = (Paap 7 3ap 7 bit ~ 3ap 7 bat + Gig — 3ap — Ly_) COS (@_ 921) (3.47) 
Tl 
Stozp (Y = (Paapm 7 3ap 7 bit — 3p 7 bi, + Gia - 3ap — Ey) 08 fo o2! + =| (3.48) 


Modulation and amplification of the data to be transmitted, is provided by the circuit 


which yields S,, (t) . 


S(t) = (Prgpm— 3ap — Lip — 3yp — Ly + Gy) C08 (@, ot + Oy (t-*)) (3.49) 
Sa (0 = Py apm COS (@_ oat + Oy (t -%)) (3.50) 

Sn = (Px in = cane — Lin + Gio) 608 [ (M92 t O.03)t + 04 (0) ] (3.51) 
Circuit at the top of the board that produces S, ,,, and S, ,,, provides an in phase and 


quadrature 10.7 MHz sinewave to the Costas loop on board 4. 


Sp = eee = eee, ) COS, Gt — oy 


T 
Sto2 (0 = (Paapm + Gi3 — 34g — Ly) ©08 (oy o2! + 5 -3| (3.52) 


The circuits in the middle of the board perform basically the same function, bi-phase 
modulation of a PN sequence onto a 21.4 MHz carrier (LO3), after some amplification and 


power splitting. 


Bet) = (Pio = Sap et cone t Or) — F-Sap — Ley) C08 (2@1530) (3.53) 
Sp (t) = Py apm 69S (OL o3t) (3.54) 

Se — CPE ee 15 _ jg COS (@ppat184)(t— %) ) G2») 
Sc(t) = Po papmy COS (@io3t + Oy (t-*)) (3.56) 
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G. BOARD 6: PN GENERATOR AND PHASE SEARCH 


Due to the lack of time, the design of board 6 is incomplete. The intention was to | 
place all the digital components on one board, thus reducing the number of boards with | 


digital ground. Digital circuits carry signals that have fast rise times and large voltages, | 


especially when compared to the weak analog signals of the RF receiver circuits. The 


concern is that a noisy digital ground may interfere with the relatively weak analog signals. | 


1. Transmitter Section 


Immediately after the message enters the board it is differentially encoded. This | 


is done to compensate for phase inversion that may occur during the detection process. | 


Naturally, the received message stream will have to be differentially decoded following 


demodulation in the receiver. The rest of the transmitter circuitry is relatively simple. | 
Once the fixed PN sequence is produced it is exclusively OR’d with the differentially | 
encoded message, and the resulting signal is delivered to board 5 for modulation and | 


transmission. The frequency synthesizer produces timing signals for synchronization of | 


the sending unit. 


2. Receiver PN Synchronization Section 


Most of the search strategies read about in References 2 and 5 conduct their search | 
for the proper PN phase using discrete shifts of the locally produced PN sequence. | 
Searches such as the one shown in Figure 3-8 may be easily programed, and handled by a | 
microprocessor. However, concerns over power consumption, and allocation of the on- : 
board processors resources lead to an implementation that would require little or no | 
processor interaction. This part of the design is original, and seeks to implement the | 


various search strategies using discrete digital components. A computer is used during the | 


evaluation phase, in order to make it easier to evaluate different search strategies. Once 
evaluated, the search strategy that provides the quickest access time could be implemented 
using a state machine to provide the parallel input to the PN sequence generators. If the 


state machine is too hard to design, then the satellites on board microprocessor should be 
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able to provide a seed sequence every 813s without taxing the microprocessors ability to 
provide other services. The search strategy shown in Figure 3-8 may be considered a 
hybrid of those shown and discussed in [Ref. 9:p. 543]. This strategy might be called a 


broken expanding window, since it is combination of the two. 


Broken / Expanding-Window 


Uncertainty 





PN Phase 






Figure 3-8 Proposed search strategy. 
A logical extension of the strategies 
shown in [Ref. 9:p. 543]. 


The goal is to produce an acquisition system that provides the shortest , 
acquisition time possible. One way to do this is to devise a search strategy that quickly 
finds the correct phase. Another method is to design a fast detection circuit, such as the one 
discussed in [Ref. 10:p. 551]. Both were to be tried, unfortunately time ran out. So the | 
faster detection circuit was not attempted. The search strategy shown in the figure above, 
is a logical extension of the strategies shown in [Ref. 10:p. 543]. This strategy takes | 
advantage of the probability that the correct phase is close to where the search is begun, and 


minimizes search time by eliminating redundancy during the search. 
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Two PN sequence generators will be used, one on line and other off line. While 
the on line PN generator is running, the off line PN generator will load 7 bits of the 127 bits 
in the sequence. When this generator is placed on line, its sequence will begin with the 7 
bits that were loaded. A phase change is introduced by shifting the 7 bits, and a search 
strategy implemented by programing both the 7 bits and the direction of the up/down 
counter. 

The buffers at the input control which PN generators will load. While the select 
lines (SAO, SA1, SBO, SB1) control whether the shift registers load, hold, orrun. The OR 
gates to the right of the PN generators are used to control which PN sequence generator is 
connected to the 8 bit shift register. The up/down counter is then used in conjunction with 
the 8 bit MUX to choose which line out of the 8 bit shift register is used. The shift register 
is clocked at twice the rate of the generated sequence, this effectively samples the sequence 
at twice the frequency. Allowing a search to be conducted in half chip steps. The 8 bit shift 
register to the right of the MUX may be clocked at either the same frequency as the first 


shift register or twice the frequency of the first shift register. If clocked at the same 


T 
frequency, the early and late sequences will be cars from punctual, and if clocked at twice 


T 
the frequency, they will be i from the generated sequence. While in search, the VCXO 


will be fixed at its fundamental frequency. Once detection has been made and tracking has 
begun, the VCXO is controlled with the feedback from the tracking loop. While tracking, 
the 7 bit string available on Port A should remain there in case the circuit loses lock and the 
search must resume. If it does loose lock, odds are that the phase is close to where it was 
when circuit lost lock. The proposed search strategy is commenced with a parallel loading 
a 7 bit seed sequence into the shift register. The 8th bit is then used to control the direction 
of the up/down counter. Looking back at the previous figure, you can see that after one 
load, the counter may count up, and after another load the counter may count down. With 


this ability, many different search strategies may be tried. 


29 


IV. CONCLUSION 


A prototype design of the complete communications subsystem was achieved during 
the course of this thesis work. This design uses a blend of analog and digital circuits to 
search and acquire the correct pseudo-noise sequence phase. The project is currently in the 
bread board phase of its design. Critical areas of the design have been tested using a crude | 
but simple RF bread board technique. Further RF bread boarding is currently under way 
using boards that have been layed-out and fabricated using a CAD/CAM system. This 
process, though very time consuming, is cleaner and should provide excellent results. | 
Unfortunately, a complete design of all the digital circuits was not accomplished. 
However, enough progress was made to start bread boarding and testing the computer 
interface and search circuitry required for system testing. 

A great deal has been learned and accomplished over the course of this design. 
However, a robust doppler compensation method and circuit remains elusive. The current 
design utilizes an orbital track prediction program, controlling a VCXO, to provide doppler 
compensation. This reliance on a prediction program should work, but is an admitted 
weakness in the design. Additional efforts to provide a more robust solution should be 


made. 
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APPENDIX A: PSEUDO NOISE SEQUENCES AND THEIR PROPERTIES 


Before going into a theoretical overview of a BPSK spread spectrum communications 
system, some fundamentals should be covered first. Most of the material on PN sequences 
was extracted from [Ref. 2:p. 1715-1718], with some re-wording for clarity. 

A. Pseudo Noise Sequences 


Pseudo noise sequences are also called maximal length sequences, or m-sequences. 


These sequences, which have a length N = 2-1, have properties that are very 
advantageous to spread spectrum communications, such as their periodic autocorrelation 


property: 
nO See Ge . | for 1<t<N-1 (A.1) 


Maximal length sequences are generated from primitive polynomials and may be 
constructed using a feedback shift register circuit. Primitive polynomials of degree ‘m’ 
are expressed in the following form 


m-— 1 


h(x) =hox "+h jx +... hx +h x +h, (A.2) 
If h(x) is primitive, it can not be factored into a product of two or more polynomials of 
lower degree. When factored down to its lowest form h, = 1 and h, = 1. Equation A.2 
will then reduce to 


-1 


h(x) =x"+h yx +. thax thx +1 (A.3) 


where ‘m’ is the span of the sequence and the number of D-flip flops used to construct the 


sequence. 
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The following example is for an m-sequence of length seven. Therefore the length of | 


the generated sequence will be N = 2'~1 = 127. The all zeros state is not a valid state. 


XOR | 





Figure A-1 Feedback shift register for a primitive 
polynomial of degree m = 7, and length 


N= gue 1 = 127, the generator polynomial 1s x + xe 


This feedback shift register generates an infinite sequence c,c,c,...c,... that satisfies the 


{ 


recurrence 


C = Care t=0,]1,... (A.4) 


t+ / 7 
In the figure above, the output is taken from the last D-flip flop in the chain. However, the — 
Same sequence may be seen, with a shift in phase, on the output of any other flip flop in | 


the chain. 


SRRBRBSL 


Figure A-2 Timing diagram for the circuit in figure 2-3 which 


was used to implement the primitive polynomial x’ +x+1. This 
timing diagram also demonstrates the shift property. 
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One way to implement this feedback shift register is as follows: 


+5V 


74HC86D 





Figure A-3 Circuit implemention of the primitive polynomial 
x’ ap ear I 
B. Properties of Maximal Length Sequences 
The following properties characterize PN sequences and provide the reader with the 
background necessary to properly apply them to an application such as spread spectrum 


communications. 


1. The Shift Property: For a PN sequence of length 2" -1, there are 2" —2 
possible shifts of the original. Each cyclic shift produces another PN sequence. The 
original and the 2" -2 cyclic shifts produces a set 5 whose size is (2° -1) x (2) -1). 

2. The Recurrence Property: Any PN sequence ce 6, satisfies the recurrence 


+...¢hc..,+c€ for <=n0sIey.. (A.5) 


tT+m = eC ra aaa ; 1~t+ 1 t 


6 


There are ‘m’ linearly independent solutions to Equation A.5, so there are ‘m’ linearly 
independent sequences in the set 6. 


3. The Window Property: If a window of width ‘m’ 1s slid along a PN sequence 


in the set 8. , each of the 2" - 1 nonzero binary m-tuples is seen exactly once. This property 


follows from the fact that the PN sequence was generated from a primitive polynomial. 


4. The Half 0’s and Half 1’s Property: A PN sequence contains 2™~" ones and 


m 


2™-'_1 zeros. Therefore there will always be one more one than there are zeros in a PN 
sequence. 

5. TheAddition Property: The sum of two sequences forms another sequence in 
the set 6. 

6. The Shift and Add Property: The sum of a PN sequence and a cyclic shift of 
it self forms another PN sequence in the set 6... 

7. The Autocorrelation Property: The autocorrelation property is the most 
important property, especially when the application is spread spectrum communications. 
This property is applied when two of the same sequences are being compared to each other. 
Such is the case with the received and locally generated sequences of a spread spectrum 


receiver. The autocorrelation function for a binary sequence c,c,c,...cy_, 18 defined by 
1 Nea . | 
Cc. Cre j 
RE a ae (A.6) 
j=0 


and for a NRZ sequence is defined by 
N-1 


l 
R.(t) = x ccC,,, ¢ = +1 (A.7) 


j=0 
The following equation may be more easily understood, and is easy to apply to both NRZ | 


and binary sequences 


Roe “ (A.8) 





where ‘A’ is the number of places where the sequences agree, ‘D’ is the number of places 
where the sequences disagree, *N’ 1s the length of the sequence, and A+D =N. 


Application of the above equations yields the following results 


] 
R.(t) = 1 when t = €N and R.(t) = A when t#EN; where € is an integer (A.9) 


Graphically the autocorrelation function looks like the following 





Figure A-4_ Plot of the autocorrelation function. The PN 
waveform is periodic with triangular pulses of width 2T, repeated 


every NT. seconds [Ref. 1:p. 407]. 


In layman's terms, the autocorrelation of two identical sequences 1s maximum when the 

two sequences are perfectly aligned in phase. If they are more than +T, apart, the 
a ] 

autocorrelation 1s a 


8. The Runs Property: In any PN sequence, 1/2 of the runs have length 1, 1/4 
have a length of 2, 1/8 have a length of 3, etc. In each case, the number of runs (repetitions) 


of zeros 1s equal to the number of runs of ones. 


a5 


APPENDIX B: LINK ANALYSIS PROGRAM AND RESULTS 


The following program and results are supplied with little explanation. Their 


inclusion is for reference only. 
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DSSS UPLINK BUDGET 
tic 172.0360 


Frx_dB= 2.0228 


Tsys = 452.0360 


Pt_Watts= 0.2000 0.2500 0.3200 0.4000 0.5000 0.6400 0.8000 
1.0000 1.2500 1.6000 2.0000 


EIRP_dBm = 35.0103 35.9794 37.0515 38.0206 38.9897 40.0618 41.0309 
42.0000 42.9691 44.0412 45.0103 


C_dBm= — -116.3669 -115.3978 -114.3257 -113.3566 -112.3875 -111.3154 -110.3463 
-109.3772 -108.4081 -107.3360 -106.3669 


N_dBm= — -108.0712 


CN_dB = -8.2957 -7.3266 -6.2545 -5.2854 -4.3163 -3.2442 -2.2751 
-1.3060 -0.3369 0.7352 1.7043 


CN_at_Det_dB = 12.7424 13.7115 14.7836 15.7527 16.7218 17.7939 
18.7630 19.7321 20.7012 21.7733 22.7424 


EbNo_dB = 15.7527 16.7218 17.7939 18.7630 19.7321 20.8042 21.7733 
22.7424 23.7115 24.7836 25.7527 


Error_Prob="0. 0 0. . 0, O00 0 Om 


PG = 21.0380 


DSSS DOWNLINK BUDGET 


rx = 171.9809 
Frx_dB= 2.0223 
Tsys = 451.9809 


Pt_Watts = 0.2000 0.2500 0.3200 0.4000 0.5000 0.6400 0.8000 
1.0000 1.2500 1.6000 2.0000 


EIRP_dBm = 23.0103 23.9794 25.0515 26.0206 26.9897 28.0618 29.0309 
30.0000 30.9691 32.0412 33.0103 


C_dBm = -116.3669 -115.3978 -114.3257 -113.3566 -112.3875 -111.3154 -110.3463 
-109.3772 -108.4081 -107.3360 -106.3669 


N_dBm = -108.0717 


CN_dB=— -8.2951 -7.3260 -6.2539 -5.2848 -4.3157 -3.2436 -2.2745 
-1.3054 -0.3363 0.7358 1.7049 


CN_at_Det_dB = 12.7429 13.7120 14.7841 15.7532 16.7223 17.7944 
i763) 19.7326 20.7017 21.7738 22.7429 


EbNo_dB = 15.7532 16.7223 17.7944 18.7635 19.7326 20.8047 
21.7738 22.7429 23.7120 24.7841 25.7532 


Error Proob= 0 0 0 0 0 0 0 0 0 0 O 


PG = 21.0380 


4] 


Jo To Go Vo Go Vo YoU Yo Go Vo Fo Vo Yo Vo Vo Go Vo Go Fo Vo Yo Vo Yo Vo Yo Fo Go Zo Fo Vo Go Vo Foo Lo Yo Vo Yo Vo Vo Yo Vo Zo Go Yo Yo Vo Yo Vo Yo Go Mo Vo Yo Vo 
% 

% 

% DSSS/BPSK UPLINK BUDGET: Amie Brown 7/14/93 

% 

Jo To To To Vo To YoU Yo Go Vo Go Vo Go Vo Vo Yo Vo Yo Vo Fo Yo To YoVo Fo Yo Yo Mo To Vo Vo Yo Yo Vo To Yo Yo Vo Yo Vo Yo Yo Yo Yo Yo Yo Mo Yo Yo Yo Yo To To Yo Vo : 
clear | 
clg 


Pt=(.2 .25 .32 .40 S50 .64 .80 1.0 1.25 1.60 2.0]; % tx pwr (watts) 
% constants 


Gain=(22.5 27 27 27]; % rx front end amplifier gains (dB) 
NF=(2 1.6 1.6 1.6]; % rx amplifier noise figures (dB) 


Gtx=12; % gain of tx antenna 

Grx=0; % gain of rx antenna 

Tant=280; % rx antenna temp due to rxvd radiation 

E=5; % worst case ground station antenna elevation in deg 
fc=436.5e6; % Carrier frequency 

r=6378.155e3; % earth radius in km 

h=480e3; % sat altitude in km 

c=2.997925e8; % speed of light 

k=1.38e-23; % Boltzmann's constant 


% converting from dB 


F=10 .A(NF/10); 

G=10 .*(Gain/10); 
Gt=10 .4(Gtx/10); 
Gr=10 .“(Grx/10); 


% calculation of the rx noise figure and noise temp 

Frx=F(1) + F(2)/G(1) + F(3)(G(1)*G(2)) + F(4)(G(1)*G(2)*G(3)); 
Trx=290*(Frx-1); 

Tsys=Tant+Trx; 


% Calculating the worst case slant range and free space loss 


d=sqrt((r+h)*2+142-2*r*(r+h)*sin(E*pi/180+asin(r/(r+h)*cos(E*pi/180)))); 
Lfs=(4*pi*fc*d/c)%2; 
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% calculating IF BW for a various data rates 


R=1:10:12000; 

R=R'; 

fchip=127*R; 

B=2*fchip; % noise bandwidth 

% calculation of the system noise temp and noise figure 

N=k*Tsys*B; 

% calculating the effective radiated pwr 

for n=1:length(Pt) 

eirp(n)=Pt(n)*Gt; 

% calculating the received signal pwr at the satellite 
Crx(n)=eirp(n)*Gr/Lfs; 


% calculating C/N and Eb/No 


CN(:,n)=Crx(n)*(ones(1:length(N)))' ./N; 
EbNo(:,n)=CN(:,n) .*B ./R; 


% calculating the probability of a bit error 


Pb(:,n)=.5*(1-erf(sqrt(EbNo(:,n))));  % fen of PSD (Eb/No) as measured at the detector input. 


end 
EN=10*log 10(EbNo); 


% calculating only for 9.84 kbps 


B1=2* 127*9.84e3; 
N1l=k*Tsys*B1; 
N2=k*Tsys*2*9.84e3; 
CN1=Crx/N1; 

CN2=Crx/N2; 
EbNo1l=CN1*B1/9.84e3; 
Pb1l=.5*(1-erf(sqrt(EbNol ))); 
PG=10*log10(127); 


% writing results to a file 


delete link 

diary link 

disp(DSSS UPLINK BUDGET) 

Trx,Frx_dB=10*log10(Frx), Tsys,Pt_Watts=Pt EIRP_dBm=10* log 1 0(eirp/le-3).... 


C_dBm=10*log10(Crx/1e-3),N_dBm=10*log10(N1/le-3),CN_dB=10*log10(CN1).... 


CN_at_Det_dB=10*log10(CN2),EbDNo_dB=10*log10(EbNo1),Error_Prob=Pb!.... 
PG 
diary off 


axis(‘square’); 

axis([0, 12000,5,35]); 

plowR,EN(, 1),R,ENCG,2),R,ENC,3),R,EN(C,4),R,ENCG,5),R,EN(,6).... 
R,EN(:,7),R,EN(:,8),R,EN(,9),R,ENC,10),R,ENG,11)),... 

tide(Bit Energy to Noise Ratio vs Information Bit Rate’)... 

xlabel( Information Bit Rate (bps)’).... 

ylabel(Eb/No (dB)’), grid, pause 

AXIS; 


plot(R,PbC,1),R,Pb(:,2),R,PbC,3),R,Pb(:,4),R,PbC,5),R,Pb(:,6).... 
R,Pb(:,7),R,Pb(:,8),R,PbC,9),R,PbC,10),R,EN(,11)).... 
tide(Probability of a Bit Error vs Information Bit Rate’)... 
xlabel(Information Bit Rate (bps)’),... 

ylabel('Pb’), grid,pause 


axis((O, 12000,0, le-6]); 
plot(R,Pb(:,1),R,PbC,2),R,PbC,3),R,PbC,4),R,PbC,5),R,Pb(:,6).... 
R,Pb(:,7),R, PbC:,8),R,PbC,9),R,Pb(:,10),R,Pb(:,11)).... 

tite( Probability of a Bit Error vs Information Bit Rate’)... 
xlabel(Information Bit Rate (bps)’).... 

ylabel(' Pb’), grid,pause 

ax1S; 


axis((O, 12000,0, le-10]); 

plot(R,Pb(:,1),R,Pb¢,2),R, PbC,3),R,PbC,4),R,Pb(-,5),R, Pb(,6).... 
R,Pb(:,7),R, PbC,8),R,PbC,9),R,PbC:,10),R,Pb(:,11)),... 

uue( Probability of a Bit Error vs Information Bit Rate’)... 
xlabel(Information Bit Rate (bps)’),... 

ylabel(’Pb’), grid,pause 

AXIS; 


axis((O, 12000,0, le-13]); 


plot(R,Pb(:,1),R,Pb¢,2),R, PbC,3),R,PbC,4),R,PbC,5),R,PbC,6),... 
R,Pb(:,7),R,PbC,8),R,Pb(:,9),R,PbC,10),R,Pb(:,11)),... 
tide(’Probability of a Bit Error vs Information Bit Rate’).... 
xlabel(Information Bit Rate (bps)’).... 
ylabel('Pb’),grid,pause 


axis; 


axis(‘Square’); 

plot(R,B,R,3e6*ones(1:length(R))).... 

tide(Null to Null Bandwidth vs Information Bit Rate’)... 
xlabel(Information Bit Rate (bps)’), 

ylabel((BW (MHz)’), grid 
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APPENDIX C: GENERALIZED IMPEDANCE CONVERTER (GIC) FILTER 
DESIGN PROGRAM 


The following program was based on notes written by Professor Micheal at the Naval | 
Postgraduate School, and was written with help from CAPT. Bing Bingham, USMC. This 
program was used to design the BPFs in Costas loop on board 4. It provides magnitude 
plots, phase plots, and component values. A component value is computed for each of the — 


variables shown in the following figure. 





Figure C-1  GIC filter topology 


The transfer function for the above topology is 


N (s) 


C.1 
D(s) (C.1) 


T(s) = 





pis) = + ( }s+ 0% (C.2) 
or : 

Each filter type has its own transfer function. A list of the filter types and its | 

corresponding variable equations and transfer function is provided in Table C-1. This | 

table is provided for clarity. It is not necessary for the proper operation of the program, | 


but the following equations are 


C= and G = (C.3) 


a | 


1 
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% GIC Realization for both Ideal and Nonideal Op Amps 

% 

% 2Jul 93 Arnie Brown 

% 

% This m-file (gic.m) uses the following “m files" writen by Bing Bingham 
% and later modified by Amie Brown: 

% 

% 1) addpoly.m adds to polynomials 

% 2) ideal.m bode plots the ideal opamp 

% 3) nonideal.m bode plots the nonideal opamp 

% 4) plotboth.m _ bode plots the nonideal over the idcal bode 
Jo 

% ALL four of these m-files are needed to run GIC.M. 

% 

% This m-file is patterned after the GIC network from page 28 of 
% Professor Micheal's notes. 


ele 
flag = 1; 
while flag == 1 


<=menu('‘Choose a Filter type’,'Low Pass’,'High Pass’,’Band Pass'’,‘Notch’,'All Pass’,"Exit GIC Program’); 


k= 
cle 
r=input(‘Enter the input resistance: '); 
fc=input(‘Enter the cutoff frequence: ‘); 
Q=input(‘Enter the quality factor (Q): '); 
c=1(r*2*pi*fc); 
yl=1/r; 
y2={c 0]; 
y3=addpoly(y2, 1/(r*Q)); 
y4=yl1; 
ys=yl; 
y6=0; 
y7=0; 
y8=y 1; 
delete LP 
diary LP 
Y 1_ohms=r, Y2_farads=c, Y 3a_farads=c, Y3b_ohms=(r*Q), Y4_ohms=r, Y5_ohms=r, Y 6=0,... 
Y7=0,Y8_ohms=r,dispCOutput is at T2’) 
diary off 
pause 
cle 
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p=menu(‘Choose IDEAL or NONIDEAL Bode plots or BOTH’, IDEAL',, NONIDEAL "BOTH ); 


cle 
if p== 
idealnum=addpoly(conv(y 1,conv(y4,y5)),addpoly(conv(y 1,conv(y5,y8)),... 
addpoly(conv(y2,conv(y3,y7)),((-1)*conv(y I ,conv(y6,y7)))))); 
idealden=addpoly(conv(y 1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.,... 
addpoly(y7,y8)))); 
ideal(idealnum, idealden,k) 


elseif p== 

wtl=input(Input the Gain-Band-Width-Product for the Op-Amps: '); 
wQ=wtl; 
nonidealnum 1l=addpoly(conv(y2,conv(y3,y7)),addpoly(conv(y 1 ,conv(y5.,... 
addpoly(y4,y8))),((-1)*(conv(y1,conv(y6,y7)))))); 

nonidealnum2=(conv(conv(y5,conv(addpoly(y 1,y3),addpoly(y4.... 
addpoly(y7,y8)))),[1 O]))/w; 

nonidealnum=addpol y(nonidealnum1] ,nonidealnum?2); 
nonidealden1=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y1,conv(y4.... 
addpoly(y5,y6)))): 
nonidealden2=((conv(conv(y 1,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2,... 
addpoly(y5,y6)))),[1 0]))/(wt1)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),{1 0]))/(wt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),(1 O]))/(wtl *wt2)); 

nonidealden=addpol y(nonidealden1,addpoly(nonidealden2.... 
addpoly(nonidealden3 ,nonidealden4))); 

nonideal(nonidealnum,nonidealden,k); 


elseif p==3 
idealnum=addpoly(conv(y1,conv(y4,y5)),addpoly(conv(y 1,conv(y5,y8)),... 
addpoly(conv(y2,conv(y3,y7)),((-1)*conv(y1,conv(y6,y7)))))); 
idealden=addpoly(conv(y 1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.... 
addpoly(y7,y8)))); 
wtl=input(Input the Gain-Band-Width-Product for the Op-Amps: '); 
wi2=wil; 
nonidealnuin1=addpoly(conv(y2,conv(y3,y7)),addpoly(conv(y 1,conv(y5.,... 
addpoly(y4,y8))),((-1)*(conv(y1,conv(y6,y7)))))); 
nonidealnum2=(conv(conv(y5,conv(addpoly(y 1,y3),addpoly(y4.... 
addpoly(y7,y8)))),[1 0)))/wt2; 
nonidealnum=addpoly(nonidealnum 1 ,nonidealnum2); 
nonidealden1=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y 1 ,conv(y4.... 
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addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y 1 ,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),{1 0]))/(wt)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),(1 0]))/(wt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),{1 0 0}))Awtl*wi2)); 

nonidealden=addpol y(nonidealden] ,addpoly(nonidealden2.... 
addpoly(nonidealden3 ,nonidealden4))); 

plotboth(idealnum, idealden,nonidealnum,nonidealden,k); 


end 


elseii kk == 2 
cle 
r=input(‘Enter the value of input resistance: ‘); 
fc=input(‘Enter the cutoff frequence: ‘); 
c=1/(r*2*pi*fc); 
Q=input(‘Enter the quality factor (Q): '); 
yl=1/r; 
y2=y1; 
y3=([c 0); 
y4=yl1; 
y5=0; 
y6=y1; 
y/=y3; 
y8=1/(r*Q); 
delete HP 
diary HP 
Y 1_ohms=r, Y2_ohmse=r, Y3_farads=c, Y4_ohms=r, Y5=0, Y6_ohms=r, Y7_farads=c,... 
Y8_ohms=(r*Q), 
dispCOutput is at T1') 
diary off 
pause 
cle 
p=menu(‘Choose IDEAL or NONIDEAL Bode plots or BOTH’, IDEAL’, NONIDEAL', BOTH); 
cle 
if p==1 
idealnum=addpoly(conv(y1,conv(y4,y5)),addpoly(conv(y3,conv(y/.... 
addpoly(y2,y6))),((-1) * conv(y3,conv(y5,y8))))); 
idealden=addpoly(conv(y 1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.,... 
addpol y(y7,y8)))); 
ideal(idealnum, idealden,k); 
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elseif p==2 
wtl=input(‘Input the Gain-Band-Width-Product for the OpAmps: '); 
wi2=witl; 
nonidealnuml=addpoly(conv(y1,conv(y4,y5)),addpol y(conv(y3,conv(y7,... 
addpoly(y2,y6))),((-1)*(conv(y3,conv(y5,y8)))))); 
nonidealnum2=(conv(conv(y7,con v(addpoly(y1,y3),addpoly(y2.... 
addpoly(y5,y6)))),{1 0)))/wul; 
nonidealnum=addpoly(nonidealnum 1 ,nonidealnum2); 
nonidealden l=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y1,conv(y4.... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y l,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),{1 O)))/(wtl)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y?2.... 
addpoly(y5,y6)))),[1 0]))/(wt2)); 
nonidealden4=((conv(conv(addpoly(y l,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),[1 0 0)))Kwtl*wi2)); 
nonidealden=addpoly(nonidealdent ,addpoly(nonidealden?.... 
addpoly(nonidealden3,nonidealden4))); 
nonideal(nonidealnum,nonidealden,k); 
elseif p== 
idealnum=addpoly(conv(y1,conv(y4,y5)),addpoly(conv(y3,conv(y7.... 
addpoly(y2,y6))),((-1) * conv(y3,conv(y5,y8))))); 
idealden=addpoly(conv(y1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.... 
addpoly(y7,y8)))); 
wtl=input(‘Input the Gain-Band-Width-Product for the OpAmps: ’); 
wQ=witl; 
nonidealnum1=addpoly(conv(y 1 ,conv(y4,y5)),addpoly(conv(y3,conv(y7,... 
addpoly(y2,y6))),((-1)*(conv(y3,conv(y5,y8)))))); 
nonidealnum2=(conv(conv(y7,conv(addpoly(y 1,y3),addpoly(y2.... 
addpoly(y5,y6)))),{1 O}))/wtl; 
nonidealnum=addpoly(nonidealnum1,nonidealnum2); 
nonidealdenl=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y 1,conv(y4.... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y l,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),[1 0]))/(wtt)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.,... 
addpoly(y5,y6)))),[1 0]))/(wt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),[(1 0 0]))/(wtl*wt2)); 
nonidealden=addpoly(nonidealden I ,addpoly(nonidealden?2.... 
addpoly (nonidealden3,nonidealden4))); 
plotboth(idealnuin,idealden,nonidealnum,nonidealden,k); 
end 
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elsen k == 3 
cle 
r=input(‘Enter the value of input resistance: '); 
fc=input(‘Enter the center frequence (Hz): '); 
B=input(‘Enter the 3dB bandwidth (Hz): '); 
c=1/(r*2* pi* fc); 
Q=fc/B; 
wl=l/r, 
y2=yl1; 
y3=[c 0); 
y4=yl, 
vo=; 
yo=yl; 
y7=1/(r*Q); 
y8=y3; 
delete BP 
diary BP 
Y 1_ohms=r, Y2_ohms=r, Y3_farads=c, Y4_ohms=r, Y 5=0, Y6_ohms=r, Y 7_ohms=(r*Q).... 
Y8_farads=c, 
disp(Output is at T1'),Q, 
diary off 
pause 
cle 
p=menu(‘Choose IDEAL or NONIDEAL Bode plots or BOTH’, IDEAL’, NONIDEAL’, ‘BOTH; 
cle 
if p== 
idealnum=addpoly(conv(y 1 ,conv(y4,y5)),addpoly(conv(y3,conv(y7.... 
addpoly(y2,y6))),((-1) * conv(y3,conv(y5,y8))))); 
idealden=addpoly(conv(y 1 ,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.... 
addpoly(y7,y8)))); 
ideal(idealnum,idealden,k); 
elseif p== 
wtl=input(Input the Gain-Band-Width-Product for the Op-Amps: '); 
wi2=wtil; 
nonidealnum1=addpoly(conv(y 1,conv(y4,y5)),addpoly(conv(y3,conv(y7.... 
addpoly(y2,y6))),((-1)*(conv(y3,conv(y5,y8)))))); 
nonidealnum2=(conv(conv(y 7,conv(addpoly(y 1,y3),addpoly(y2.... 
addpoly(y5,y6)))),[1 0)))/wtl; 
nonidealnum=addpoly(nonidealnum 1 ,nonidealnum?2); 
nonidealden 1 =addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y 1 ,conv(y4.... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y1,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 


addpoly(y5,y6)))),[1 O))/(wu)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),{1 0)))/(wt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),{1 0 0)))/(wtl*wt2)); 
nonidealden=addpoly(nonidealden!l ,addpoly(nonidealden2.... 
addpoly(nonidealden3 nonidealden4))); 
nonideal(nonidealnum,nonidealden,k); 
elseif p== 
idealnum=addpoly(conv(y 1,conv(y4,y5)),addpoly(conv(y3,conv(y7.... 
addpoly(y2,y6))),((-1) * conv(y3,conv(y5,y8))))); 
idealden=addpoly(conv(y1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.... 
addpoly(y7,y8)))); 
wtl=input(Input the Gain-Band-Width-Product for the OpAmps: '); 
wQ=wtl; 
nonidealnum 1l=addpoly(conv(y 1,conv(y4,y5)),addpoly(conv(y3,conv(y7,... 
addpoly(y2,y6))),((-1)*(conv(y3,conv(y5,y8)))))); 
nonidealnum2=(conv(conv(y7,conv(addpoly(y 1,y3),addpoly(y2.... 
addpoly(y5,y6)))),[1 0)))/wal; 
nonidealnum=addpoly(nonidealnum!] ,nonidealnum2); 
nonidealden 1=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y 1,conv(y4.... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y 1,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),(1 0)))/(wtl)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),[1 O)))(wt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),(1 0 0J))/(wul *wi2)); 
nonidealden=addpoly(nonidealden! ,addpoly(nonidealden2,... 
addpoly(nonidealden3,nonidealden4))); 
plotboth(idealnum, idealden,nonidealnum,nonidealden,k); 


end 


elseif k == 
cle 
r=input(‘Enter the value of input resistance: '); 
fc=input(‘Enter the center frequence: '); 
B=input(Enter the 3dB bandwidth: '); 
c= 1/(r*2*pi* fc); 
Q=fc/B; 
yl=1/r; 
y2=y1; 
y3=[c 0]; 


y4=y 1; 
y5S=y1; 
y6=0; 
WENGE 
y8=1/(r*Q); 
delete N 
diary N 
Y 1_ohms=r, Y2_ohmsz=r, Y3_farads=c, Y4_ohms=r, Y5S_ohms=r, Y6=0, Y7_farads=c.... 
Y8_ohms=(r*Q), 
dispCOutput is at T2'),Q 
diary off 
pause 
CIC 
p=menu('Choose IDEAL or NONIDEAL Bode plots or BOTH’, TIDEAL',, NONIDEAL', BOTH’); 
cle 
ip==1 
idealnum=addpoly(conv(y 1 ,conv(y4,y5)),addpoly(conv(y 1,conv(y5,y8)),... 
addpoly(conv(y2,conv(y3,y7)),((-1)*conv(y 1,conv(y6,y7)))))); 
idealden=addpoly(conv(y 1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.... 
addpoly(y7,y8)))); 
ideal(idealnum,idealden,k); 
elseif p== 
wtl=input(‘Input the Gain-Band-Width-Product for the OpAmps: '); 
w=wtl; 
nonidealnuin1=addpoly(conv(y2,conv(y3,y7)),addpol y(conv(y 1 ,conv(y5.... 
addpoly(y4,y8))),((-1)*(conv(y1l,conv(y6,y7)))))); 
nonidealnum2=(conv(conv(y5,conv(addpoly(y 1,y3),addpoly(y4.... 
addpoly(y7,y8)))),[1 0)))/wt2; 
nonidealnum=addpol y(nonidealnum 1 ,nonidealnum2); 
nonidealden l=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y l,conv(y4.,... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y 1 ,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),[1 0))wt1)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),[1 0})Awt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),[1 0)))/(wtl *wt2)); 
nonidealden=adupol y(nonidealden1 ,addpoly(nonidealden?.... 
addpol y(nonidealden3 nonidealden4))); 
nonideal(nonidealnum, nonidealden,k); 
elseif p==3 
idealnum=addpoly(conv(y 1 ,conv(y4,y5)),addpoly(conv(y 1,conv(y5,y8)).... 
addpoly(conv(y2,conv(y3,y7)),((-1)*conv(y 1,conv(y6,y7)))))); 
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idealden=addpoly(conv(y 1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.,... 
addpoly(y7,y8)))); , 
wtl=input(Input the Gain-Band-Width-Product for the OpAmps: '); 
wQ=wil; 
nonidealnum1=addpoly(conv(y2,conv(y3,y7)),addpoly(conv(y 1,conv(y5.... 
addpoly(y4,y8))),((-1)*(conv(y1,conv(y6,y7)))))); 
nonidealnum2=(conv(conv(y5,conv(addpoly(y1,y3),addpoly(y4.... 
addpoly(y7,y8)))),(1 0)))/wt2; 
nonidealnum=addpoly(nonidealnum 1] ,nonidealnum2); 
nonidealden l=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y1,conv(y4.... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y 1,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),[1 0]))/(wt1)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),{1 0]))/(wt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),[1 O}))/(wtl *wt2)); 
nonidealden=addpoly(nonidealden] ,addpoly(nonidealden2.... 
addpoly(nonidealden3 ,nonidealden4))); 
plotboth(idealnum, idealden,nonidealnum,nonidealden,k); 
end 


elseif k =——) 
cle 
r=input(Enter the value of input resistance: '); 
fc=input(‘Enter the center frequence: '); 
Q=input(‘Enter the value of Quality factor (Q): '); 
cal /(rsZ piste): 
yl=1/r; 
y2=yl; 
y3=[c 0); 
y4=yl; 
yo=y 1; 
y6=0; 
y7=y3; 
y8=1/(r*Q); 
delete AP 
diary AP 
Y 1_ohms=r, Y2_ohms=r, Y3_farads=c, Y4_ohms=r, Y5_ohms=r, Y6=0, Y7_farads=c.,... 
Y8_ohins=1/y8, 
disp(Output is at T1') 
diary off 
pause 
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cle 
p=menu('Choose IDEAL or NONIDEAL Bode plots or BOTH ', IDEAL',,NONIDEAL', 'BOTI'); 
cle 
if p==1 
idealnum=addpoly(conv(y 1,conv(y4,y5)),addpoly(conv(y3,conv(y7.... 
addpoly(y2,y6))),((-1) * conv(y3,conv(y5,y8))))); 
idealden=addpoly(conv(y1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.... 
addpoly(y7,y8)))); 
ideal(idealnum, idealden,k); 
elseif p==2 
wtl=input(‘Input the Gain-Band-Width-Product for the Op-Amps: '); 
wQ=wtl; 
nonidealnum 1=addpoly(conv(y 1 ,conv(y4,y5)),addpoly(conv(y3,conv(y 7... 
addpoly(y2,y6))),((-1)*(conv(y3,conv(y5,y8)))))): 
nonidealnum2=(conv(conv(y7,conv(addpoly(y 1,y3),addpoly(y2.... 
addpoly(y5,y6)))),{1 0]))/wel; 
nonidealnum=addpoly(nonidealnum1,nonidealnum2); 
nonidealden l=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y 1,conv(y4.... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y 1 ,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),(1 0)))/(wt1)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),[1 0]))/(wt2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)),... 
addpoly(y2,addpoly(y5,y6)))),[1 0 0]))(wtl*wt2)); 
nonidealden=addpoly(nonidealden 1 ,addpoly(nonidealdenz2.... 
addpoly(nonidealden3 ,nonidealden4 ))); 
nonideal(nonidealnum,nonidealden,k); 
elseif p== 
idealnum=addpoly(conv(y 1,conv(y4,y5)),addpoly(conv(y3,conv(y 7... 
addpoly(y2,y6))),((-1) * conv(y3,conv(y5,y8))))); 
idealden=addpoly(conv(y1,conv(y4,addpoly(y5,y6))),conv(y2,conv(y3.... 
addpoly(y7,y8)))); 
wtl=input(‘Input the Gain-Band-Width-Product for the Op-Amps: '); 
wQ=wil; 
nonidealnum1=addpoly(conv(y1,conv(y4,y5)),addpoly(conv(y3,conv(y 7... 
addpoly(y2,y6))),((-1)*(conv(y3,conv(y5,y8)))))); 
nonidealnum2=(conv(conv(y7,conv(addpoly(y 1,y3),addpoly(y2.... 
addpoly(y5,y6)))),[1 0]))/wtl; 
nonidealnum=addpoly(nonidealnum 1 nonidealnum2); 
nonidealden1=addpoly(conv(y2,conv(y3,addpoly(y7,y8))),conv(y 1 ,conv(y4.,... 
addpoly(y5,y6)))); 
nonidealden2=((conv(conv(y1,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
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addpoly(y5,y6)))),{1 0}))/(wt1)); 
nonidealden3=((conv(conv(y3,conv(addpoly(y4,addpoly(y7,y8)),addpoly(y2.... 
addpoly(y5,y6)))),(1 0)))/(wi2)); 
nonidealden4=((conv(conv(addpoly(y 1,y3),conv(addpoly(y4,addpoly(y7,y8)).... 
addpoly(y2,addpoly(y5,y6)))),(1 0 0)))/(wtl *wt2)); 
nonidealden=addpoly(nonidealden1 ,addpoly(nonidealden?.... 
addpoly(nonidealden3 nonidealden4))); 
plotboth(idealnum,idealden,nonidealnum,nonidealden,k); 


end 
elseif k == 6; 
break; 
elseif k == 99; 
keyboard; 
end 


end 


funcuon y=addpoly(a,b) 

Jo 

% EC4100 Bing Bingham 

Jo 

% Network Theory Dr. Sherif Michael 

Jo 

% 21 Aug 92 

Jo 

% This function adds two polynomials together after ensuring 
% that they have same length so as not to lose their correct 
% exponent value. Different length polynomials are zero 
% padded. 

% 

% 


if length(a) >= Iength(b) 
b=[zeros(1,length(a) - length(b)) b]; 
else 
a=(zeros(1,length(b) - length(a)) al]; 


end 


y = atb; 
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function y=ideal(idealnum,idealden,k) 
J 
J 2 Jul 93 Amie Brown 
Jo 
% This function draws the bode plots for an ideal op-amp. 
Jo 
% “k" is obtained from the opening menu and is passed in. 
% 
cle 
fl=input(‘Enter the BEGINNING frequency for the bode plot in powers of 10 (1Hz=0): '); 
f2=input(‘Enter the ENDING frequency for the bode plot in powers of 10 ((QOMHz=8): '); 
f=logspace(fl,f2,500); 
cle;clg 
[idealmag,idealphase]=bode(idealnum ,idealden,2* pi*f); 
idealmag=idealmag/max(idealmag); 
if k== 
subplot(211) 
x=find(idealmag>. 707); 
semilogx(f,20*log10(idealmag),f(1:max(x)),-3 *ones(1:max(x)).... 
[f(max(x)) f(max(x))J,[min(20*log10(idealmag)) -3],'g'), grid; 
udle(Magnitude Plot fora LOW PASS GIC Filter using IDEAL Op-Amps)); 
elseif k==2 
subplot(211) 
x=find(idealmag>.707); 
semilogx(f,20*log10(idealmag).... 
f(x),-3 *ones(x),... 
[f(min(x)) f(min(x))],[min(20*log10(idealmag)) -3],'g'), grid; 
uidle(Magnitude Plot fora LOW PASS GIC Filter using IDEAL Op-Amps)); 
elseif k== 
subplot(211) 
x=find(idealmag>.707); 
semilogx(f,20*log10(idealmag).... 
[f(min(x)) f(min(x))],[min(20*log 10(Gdealmag)) -3}.,... 
[f(max(x)) f(max(x))J,[min(20*log10(idealmag)) -3],'g’,... 
f(x),-3*ones(1:length(x)),'g'), grid; 
ude(‘Magnitude Plot fora BAND PASS GIC Filter using IDEAL Op-Amps’); 
elseit K== 
subplot(211) 
semilogx(f,20*log10(idealmag)), grid; 
uue(Magnitude Plot fora NOTCH GIC Filter using IDEAL Op-Amps’); 
elsetl k== 
subplot(211) 
semilogx(f,20*log10(idealmag)), grid; 


ay 


uue(‘Magnitude Plot for a ALL PASS GIC Filter using IDEAL Op-Amps’); 
end 
xlabel(‘Frequency (Hz)');ylabel(‘Gain (dB)’); 
subplot(212); 
semilogx(f,idealphase), grid; 
if k== 

title(Phase Plot fora LOW PASS GIC Filter using IDEAL Op-Amps’); 
elseif k==2 

tide(‘Phase Plot for a HIGH PASS GIC Filter using IDEAL Op-Amps’); 
elseif k==3 

uule(‘Phase Plot fora BAND PASS GIC Filter using IDEAL Op-Amps'’); 
elseif k== 

ule(‘Phase Plot fora NOTCH GIC Filter using IDEAL Op-Amps’); 
elseif k==5 

ude(‘Phase Plot for a ALL PASS GIC Filter using IDEAL Op-Amps’); 
end 
xlabel(Frequency (Hz)’); 
ylabel('Phase (degrees)’); 
pause; 


function y=nonideal(nonidealnum,nonidealden,k) 
% 
% This function draws the bode plots for an nonideal opamp 
% 
% "k" is obtained from the opening menu and is passed in. 
% 
clc 
fl=input(‘Enter the BEGINNING frequency for the bode plot in powers of 10 (1Hz=0) '); 
f2=input(‘Enter the ENDING frequency for the bode plot in powers of 10 (IOOMHz=8) '); 
f=logspace(fl ,f2,500); 
cle;clg 
[nonidealmag,nonidealphase]=bode(nonidealnum,nonidealden,2* pi*f); 
nonidealmag=nonidealmag/max(nonidealmag); 
if k== 
subplot(211) 
x=find(nonidealmag>.707); 
semilogx(f,20*log 1 O(nonidealmag),f(1:max(x)),-3*ones(1:max(x)),... 
[f(max(x)) f(max(x))],[min(20*log10(nonidealmag)) -3),'g'),grid; 
tile(Magnitude Plot fora LOW PASS GIC Filter using NONIDEAL Op-Amps)); 
elseif k== 
subplot(211) 
x=find(nonidealmag>.707); 
semilogx(f,20*log10(monidealmag).... 


f(min(x):length(nonidealimag)),-3 *ones(min(x):length(nonidealimag)).... 

[f(min(x)) f(min(x))], (min(20*log 10(nonidealmag)) -3)],'g'), grid; 

title’Magnitude Plot fora LOW PASS GIC Filter using NONIDEAL Op-Amps'); 
elseif k==3 

subplot(211) 

x=find(nonidealmag>.707); 

semilogx(f,20* log 10(nonidealmag).... 

[f(min(x)) f(min(x))],[(min(20*log 10(nonidealmag)) -3}.... 

[f(max(x)) f(max(x))],{(min(20*log1 O(nonidealmag)) -3],'g’,... 

f(x),-3*ones(1:length(x)),'g’), grid; 

tide(Magnitude Plot fora BAND PASS GIC Filter using NONIDEAL Op-Amps’); 
elseif k==4 

subplot(211) 

semilogx(f,20*log 10(nonidealimag)), grid; 

ude(Magnitude Plot for a NOTCH GIC Filter using NONIDEAL Op-Amps’); 
elseif k==5 

subplot(211) 

semilogx(f,20*log10(nonidealmag)), grid; 

utle?Magnitude Plot for a ALL PASS GIC Filter using NONIDEAL Op-Amps’); 
end 
xlabel(‘Frequency (Hz)’);ylabel(‘Gain (dB)’); 
subplot(212); 
semilogx(f,nonideal phase), grid; 
if k==1 

tide(Phase Plot for a LOW PASS GIC Filter using NONIDEAL Op-Amps’); 
elseif k==2 

tite(Phase Plot fora HIGH PASS GIC Filter using NONIDEAL Op-Amps’); 
elseif k==3 

title(Phase Plot fora BAND PASS GIC Filter using NONIDEAL Op-Amps’; 
elseif k==4 

title(’Phase Plot for a NOTCH GIC Filter using NONIDEAL Op-Amps’); 
elseif k==$ 
end 
xlabel(Frequency (Hz)’); 
ylabel(‘Phase (degrees)’); 
pause; 


funcuon y=plotboth(num,DEN,NUM,DENN,k) 

% 

% 2 Jul 93 Amie Brown 

% 

% This function draws the bode plots for an both ideal and non-ideal op-amps. 
Jo 
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% “k" is obtained from the opening menu and is passed in. 

Jo 

cle 

f1=input(‘Enter the BEGINNING frequency for the bode plot in powers of 10 (1Hz=0) '); 
f2=input(‘Enter the ENDING frequency for the bode plot in powers of 10 (IQOMHz=8) °); 
f=logspace(f1,f2,500); 

cle 

clg 

[idealmag,idealphase]=bode(num,DEN,2* pi*f); 

idealmag=idealmag/max(idealmag); 

[nonidealmag ,nonidealphase]=bode(NUM,DENN,2* pi* f); 
nonidealmag=nonidealmag/max (nonidealmag); 


if k==1 

x=find(idealmag>.707); 

subplot(211); 

semilogx(f,20*Ilog 1 0(idealmag),f,20*log 10(monidealmag),'b’,... 

[f(max(x)) f(max(x))],[min(20*log10(monidealmag)) -3],... 

f(1:max(x)),-3*ones(1:max(x)),'g'), grid; 

tide(LOW PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps)); 
elseitk=-=2 

x=find(idealmag>.707); 

subplot(211); 

semilogx(f,20*log 1 0(idealmag),f,20* log 10(nonidealmag), b’.... 

f(min(x):length(nonidealmag)),-3*ones(min(x):length(nonidealmag)).... 

[f(min(x)) f(min(x))],[min(20*log 10(nonidealmag)) -3],'g'), grid; 

ude( HIGH PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps’); 
elseif k== 

x=find(idealmag>.707); 

subplot(211); 

semilogx(f,20*log 1 0(idealmag),f,20* log 10(nonidealmag), b’,... 

[f(min(x)) f(min(x))],[min(20*log 10(nonidealmag)) -3],... 

[f(max(x)) f(max(x))], [min(20* log 10(nonidealmag)) -3],'g’,... 

f(x),-3*ones(1:length(x)),'g'), grid; 

tile(BAND PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps’); 
elseif k==4 

x=find(idealmag>.707); 

subplot(211); 

semilogx(f,20* log 10(idealinag),f,20* log 10(nonidealinag),'b’), grid; 

title(‘4th Order BAND PASS GIC Filter using NON-IDEAL Op-Amps’); 

xlabel(Frequency (Hz)'); 

ylabel('Gain (dB)’); 

subplot(212) 
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semilogx(f,idealphase,f,nonideal phase,'b'), grid, pause,clg; 
nonidnum=conv(num,NUM);nonidden=conv(DEN,DENN); 
[nonidmag,nonidphase]=bode(nonidnum, nonidden, 2* pi*f); 
nonidmag=nonidmag/max(nonidmag); 
semilogx(f,20*log10(nonidmag).... 
[f(min(x)) f(min(x))],[min(20*log 10(nonidmag)) -3],... 
[f(max(x)) f(max(x))],[min(20*log10(nonidmag)) -3],'g’,... 
f(x),-3 *ones(1:length(x)),'g'), grid; 
subplot(212) 
semilogx(f,nonidphase), grid,pause,clg; 
elseif k==° 
subplot(211) 
semilogx (f,20*log10(idealmag),f,20*log 1 0(nonidealmag),’b'’) 
ude(NOTCH GIC Filter using both IDEAL & NON-IDEAL Op-Amps’); 
elseif k== 
subplot(211) 
semilogx(f,20*log 10(idealmag),f,20*log 10(nonidealmag),'b') 
ude(ALL PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps'); 
end 
xlabel(Frequency (Hz)’); 
ylabel(‘'Gain (dB)’'); 
subplot(212) 
semilogx(f,idealphase, f,nonidealphase,'b’), grid; 
if k==] 
ude(LOW PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps’); 
elseif k==2 
uudeC HIGH PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps’); 
elseif k== 
ude(BAND PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps’); 
elseif k==4 
ude(NOTCH GIC Filter using both IDEAL & NON-IDEAL Op-Amps'’); 
elseif k==5 
uide('ALL PASS GIC Filter using both IDEAL & NON-IDEAL Op-Amps’); 
end 
xlabel(Frequency (Hz)'); 
ylabel(‘Phase (degrees)'); 
pause; 
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APPENDIX D: DESCRIPTION OF CIRCUIT SCHEMATICS 
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Figure D-1_ LPKF Circuit Board Plotting (CBP) 
CAM System. 


This machine mills around pads and traces to isolate them on a copper clad board. The; 
excess copper of both the component and solder sides have been used to form a signal’ 
ground plane. 
A. BOARD 1: TRANSMITTER AND RECEIVER FRONT END 

Board 1 was not milled because all the components have SMA connectors on their 
input and output. Connections made in this way avoided high frequency circuit board 


design problems, and made component placement less critical. Though the design is fairly! 


straight forward, there are still a few areas of concern such as the PIN diode switch (U2), 


preselection filter (U3), and the power amplifier (U4). 


== 





Figure D-2 RF bread board of the transmitter and 
receiver front end. 


The PIN diode switch (U2) allows the user to select between two separate transmitters 
and two separate receivers. However, the component chosen has some limitations. The 
first limitation is a 1 dB compression point of 19 dBm and a maximum RF power of 28 
dBm. Above 19 dBm the diodes will start to distort the signal, and above 28 dBm damage 
to the device may occur. Second limitation, it is unclear from the data sheet whether the 
Switch is absorptive or reflective. If the device is absorptive, a port that is switched off will 
appear as a 50Q load, and if it is reflective it may appear as a short [Ref. 8:p. 13-3 to 13-4]. 
Obviously, the absorptive case is the desirable one in this application. Of course, this 
should be measured before using the switch in a live circuit. The next area of concern is 
the high power amplifier (U4). This device was chosen because it was readily available 
from one of the labs here, and its maximum output power is not as likely to damage the PIN 
diode switch (U2). The last area of concern is the preselection filter (U3). This filter is 


intended to limit the band of frequencies that are allowed to enter the low noise amplifier. 
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It may not be necessary, but was included to eliminate strong out of band carriers that could 
possibly saturate the low noise amplifier. 


B. BOARD 2: RECEIVER IF AMPLIFICATION, AGC, AND DETECTOR 





Figure D-3_ RF bread board of the 
Rx IF amplification, AGC, and 
detection circuits of board 2. 







1. Rx IFAmplification and AGC Circuit 

The purpose of this circuit is to provide a relatively constant amplitude | 

wideband signal to the detector, tracker, and demodulator circuits. U1, U2, U4, and U9 | 
supply the gain, while U6 attempts to maintain a constant amplitude signal. U6 is an AGC | 


amplifier with a O-5 volt control voltage. Its feedback control voltage is developed with 


U11, U12, and U7. U11 provides a portion of the signal power to U12, the analog level © 





detector, without causing reflections on the main signal path. U12 contains both reference 


and detector diodes, as shown in the figure below. 





RFw 


Reference 


c------ - 
uu 
o 


Figure D-4 Schematic diagram of the analog 
level detector (U12) [Ref. 6:p. (5-10)]. 


These diodes are provided a bias current to increase their sensitivity. The biasing is 
cancelled out by summing the detectors output with the reference. The AGC amplifier is 
capable of providing both gain and attenuation, leading to a range of over 36 dB. When 
the signal power applied to the input of the AGC amp is 2.5 dBm, the feedback path 
should provide 5 volts back to the AGC control pin. When this occurs, the AGC amplifier 
will attenuate the signal by 13 dB. 
2. Detector Circuit 

The purpose of this circuit is to detect when the received and locally produced 

PN sequences come within half a chip of being in phase. Figure D-5 is a schematic diagram 


of the threshold detector (U3) used on board 2. When coarse phase alignment is achieved, 
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the signal will de-spread from a wideband signal to a narrowband signal. Once this occurs, 


the signal power applied to the input of U3 should be sufficient to trigger a detection. 


Vee (+8 to +20 VOC 


THRESHOLD ADJUST 
AT +15 VOC BIAS 


SENSITIVITY 
dBm (typ) 








Figure D-S Schematic diagram of the threshold 

detector (U3) [Ref. 6:p. (5-16)]. 
Unlike the analog detector this detector does not have to be externally biased, its 
sensitivity 1s determined by the value of the resistance seen by the threshold adjust pin. 
| 


should be kept as narrow as possible. However, as the bandwidth goes down, the required ! 


To reduce the possibility of a false detection due to noise, the bandwidth of the BPF 


Q for the filter goes up. This is why a crystal filter was used. Unfortunately, the crystal 


filter and its impedance matching circuitry exhibited a higher insertion loss than originally | 


anticipated. The additional insertion loss of the crystal BPF was then compensated for 


\ 


with U1L3. ! 





Each crystal filter is a two pole BPF with an input and output impedance of 3kQ | 


and 2pf. To achieve a steeper role off, two crystal filters are cascaded together to create a 


four pole filter. The input and output must be impedance matched to the source and load, 
each of which is 5082. Since the signal is only 15kHz wide at this point, a narrowband | 
impedance matching circuit is all that is necessary. A narrowband impedance matching | 


circuit has fewer components than the wideband impedance matching circuit and is easier 


to build. Initial values for the impedance matching circuit were computed using a free CAD 


program from Hewlett Packard called APPCAD. 


Hewlett-Packard 


JjXi --—o 
2 pf 
Ri N 3kn 
/ 
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ki Xi 
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oO jxb [-—o o—{ jXo a 
15.6 pH 00 
/ 
JXa_ 12.6 pr soQ \ Ro 
/ 
) Matching Ckt Load 
177 107 
Xa Xb Ro Xo 
-1177.67 1034.10 39 .60 9.00 
925.57 -1034.09 Tmag Pang 
6.06 30.0 
Zo = 50 


Impedance Mode 


(F2)=Compute [F3]=Mode [ESCI=Quit 


Figure D-6 Snap shot of the monitor screen 
while using APPCAD to design the impedance 
matching circuits. 


These values were used as a Starting point during a crude RF breadboard of the circuit. A 


variable inductor and variable capacitor were used to determine the best values. As it 


turned out, the inductor had less of an effect on the circuit than the capacitor. Soa 


variable capacitor will be used in the final bread board to minimize insertion loss of the 


impedance matching circuit. 
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C. BOARD 3: PNSEQUENCE TRACKING 





Figure D-7 Partially completed RF bread 
board of the PN sequence tracking circuit 
on board 3. 





The purpose of this circuit is to provide an error voltage to the VCXO of the locally | 
produced PN sequence generator to allow the locally produced PN sequence to track the | 
phase (or epic) of the received PN sequence. In this design, this is accomplished with an | 
implementation of the non-coherent double dither loop (DDL) [Ref. 5:p. 188-189]. 

Turning to the schematic of board 3, you may notice a few of the same components : 


and circuits used in the AGC and detector circuit diagrams. This is because each arm of | 
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the analog detectors is expected to be - 6 dBm yielding an output of approximately 200mv. 


This estimate of the output 1s based on figure D-8. 


Detected Ouiput 
ve. Power Input Ree 





-20 -15 -10 - 0 «5 
Power inpui, dBm 


Figure D-8 Graph showing the detected 

output of the analog detectors vs input 

power [Ref. 6:p. (5-11)]. 
The output from each detector, after the reference and detected voltages are summed 
(U17, U18), then subtracted (U20, U21), to develope and error voltage. The error 
voltage is amplified (U16) and low pass filtered before it is used to control the 5 MHz 
HCMOS VCXO on board 6. An indication of lock is developed using a window 
comparator (U15) which senses when the loop feedback voltage’s magnitude exceeds that 
measured when the loop is tracking. An inverter should be used on board 6 to make this 
signal a positive voltage when the loop is locked and tracking. 
D. BOARD 4: CARRIER TRACKING LOOP AND DEMODULATOR 

The purpose of this circuit is to align the phase of the locally produced carrier with 

that of the received suppressed carrier by producing an error voltage that controls the 10.7 
MHz sinewave VCXO on board 5. This is necessary to properly demodulate a BPSK 
signal, which may only be demodulated coherently. The circuit is an implementation of a 
Costas loop. This board has not been fabricated because a decision was made to include 
demodulator circuitry, and that circuitry still remains to be designed. As you can see, some 
of the circuitry is the same as that used in the previous schematics, and will not be repeated 


here. 
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The LPFs that appear in each arm (U9, U10, U12, U13) were implemented using the 
GIC topology as discussed in Appendix C. The output of each LPF when multiplied 


together (U14) yields an error signal. 


(Xy~— Xz ) (¥y~Y¥2) (Bus Re) 
oe 10V mi )*> 
1kNs A1,A2< 100kN 





Figure D-9 Four-quadrant analog multiplier 
used as a phase detector on board 4 [Ref. 13:p. 
(2-47)]. 


This error signal represents the phase difference between the two signals, because the two 
signals are the same frequency and in quadrature. The remaining parts of the circuit such 
as the amplifier, loop filter, and window comparator are identical to those used on board 3 


and discussed in the previous section. 
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E. BOARD 5: PN PHASE MODULATION AND AMPLIFICATION 





Figure D-10 RF bread board of the PN phase 
modulation and amplification circuitry on board 5. 


The purpose of this circuit 1s to prepare binary data and PN sequences for mixing. 
This is accomplished by first bi-phase modulating them on to a carrier, then amplifying the 
resulting BPSK signal to the level required. Because a number of circuits required bi-phase 
modulation, it was convenient to place all of them on the same board. 

Before arriving on this board the binary data and PN sequences are converted from 
voltage to current, because the bi-phase modulators (U7, U17, U18, U19, U20) are driven 


with + 10 milliampere. This may be done several ways, one of which 1s shown below 





Figure D-11 Circuit used to convert a TTL 
signal to the +10ma drive current required by 
the bi-phase modulators. 


16: 


To function properly, the carrier power applied to the input of the bi-phase modulators 
should not exceed -3 dBm. Once functioning, the resulting BPSK signal is then amplified | 


to the level required for mixing on the boards that follow. : 


+15V 


Output 
Input; 





AGND AGND -15V 


Figure D-12 Typical non-inverting amplifier used on board 5 to 
amplify a signal to the level required for mixing. 


U1 of figure D-9 is a current feedback operational amplifier. A current feedback op- 
amp was selected because its bandwidth is nearly independent of the closed loop gain.. 
Unlike a voltage feedback amplifiers, whose bandwidth drops rapidly as closed loop gain| 


is increased. This op-amp is very sensitive to the value of R1 as shown in figure D-9. 


Table D-1: Feedback resistor values recommended by 
Harris Semiconductor to maintain both bandwidth and stability 
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Initial attempts to build this circuit on a crude RF bread board failed, because the value of 
R, was too high. After several calls to Harris Semiconductor, the engineer who designed 
the chip provided the information shown in Table D-1. The amplifier worked great once 
the proper value of R,. was installed. The fact that a current feedback op-amp Is sensitive 
to the value of the feedback resistor was not obvious at the time, because a voltage 


feedback amplifier is what is generally discussed in text books. 


F. BOARD 6: PN GENERATION AND PHASE SEARCH CIRCUITS 





Figure D-13 Digital bread board of the PN phase 
search circuit shown in the schematic diagram of 
Board 6. 

The purpose of this circuit 1s to conduct a search for the proper phase of the locally 
produced PN sequence. Many search strategies have been written about as discussed in 
chapter three. This design attempts to implement and allow the evaluation of these 
different search strategies, by allowing the user to program the search strategy to be run. 
Once the evaluation phase is complete, the desired search strategy may be implemented 
with a state machine. The state machine would take the place of the computer in the 
evaluation phase. 

This design is completely experimental, and is not based on anything seen or read. 


Fortunately, because it is digital, it can be easily adapted without having to produce another 


board, as could occur with the RF bread boards. A detailed discussion of the actual circuit | 


may be fruit less at this point, because the schematic 1s incomplete. 
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Figure D-14 Layout diagram for Board 2 the Rx IF 
Amplification, AGC, and Detection circuits. 
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Figure D-15 Component side of Board 2. 








Figure D-16 Solder side of Board 2. 
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Figure D-17 Layout diagram for Board 3 the PN Sequence 
and Tracking circuits. 
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Figure D-18 Component side of Board 3. 
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Figure D-19 Solder side of Board 3. 
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Figure D-20 Layout diagram for Board 4 the Carrier 
Tracking Loop and Demodulator circuits. 


89 








o- oa 





230 
0240 
250 
260 
2790 


on fm OM & © NY = 
Oo 
o 





eooooooooo0o0o00c0>d: 





15 


as) 











08 120 160 


= 
& 


Layout diagram for Board 5 the Modulati 
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Figure D-22 Component side of Board 5. 
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Figure D-23 Solder side of Board 5. 
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